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Abstract. We have studied the relationships between the radio continuum (RC) and CO emission for a set of 
galaxies selected from the BIMA Survey of Nearby Galaxies. We find that the global CO-RC correlation is as 
tight as the global FIR-RC correlation for the 24 galaxies studied. Within 9 galaxies with ^6” CO and RC data 
available, the CO and RC emission is as tightly correlated as its global value; the radially averaged correlation 
is nearly linear, extends over four order of magnitude and holds down to the smallest linear resolution of the 
observations, which is ^ 100 pc. We define qco/nc the log of the ratio of the CO to RC flux as a way to 
characterize the CO-RC correlation. Combining 6" pixel-by-pixel comparisons across all sources yields an average 
small-scale correlation of ?co/rc = 1-1 i 0.28; that is, the spatially resolved correlation has a dispersion that 
is less than a factor of 2. There are however systematic variations in the CO/RC ratio; the strongest organized 
structures in geo /rc tend to be found along spiral arms and on size scales much larger than the resolution of the 
observations. We do not measure any systematic trend in CO/RC ratio as a function of radius in galaxies. The 
constancy of the CO/RC ratio stands in contrast to the previously measured decrease in the FIR/RC ratio as a 
function of radius in galaxies. We suggest that the excellent correlation between the CO, RC and FIR emission 
in galaxies is a consequence of regulation by hydrostatic pressure; this model links all three emissions without 
invoking an explicit dependence on a star formation scenario. 
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1. Introduction 

The correlation between the global far-infrared (FIR) and 
centimeter-wavelength radio continuum (RC) emissions in 
galaxies is at the same time one of the most robust and 
one of the most puzzling relations in extragalactic re¬ 
search. On global scales, the FIR and RC emissions are 
linearly correlated over five orders of magnitude of lumi¬ 
nosity, with an RMS scatter of less than a factor of 2 
(Condon 1992; Yun, Reddy, & Condon 2001). The rela¬ 
tion holds across a remarkable variety of galaxies to in¬ 
clude dwarf galaxies (Leroy et al. 2004), normal barred 
and unbarred spiral galaxies of all Hubble types, irreg¬ 
ular galaxies, E/SO galaxies with ongoing star formation, 
starbursts, Seyferts and radio-quiet quasars (Condon 1992 
and references therein), out to the most extreme starform¬ 
ing galaxies in the early Universe, at z = 1 and beyond 
(Appleton et al. 2004). 
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Although both the FIR and RC emission are thought 
to originate as a result of star formation in molecular 
clouds, the emission mechanisms for the two processes are 
completely different. The FIR arises from dust heated by 
newly-born massive stars, still embedded in their parent 
envelopes. The RC is, for the most part, nonthermal syn¬ 
chrotron emission that arises from the interaction of rela¬ 
tivistic electrons with the ambient magnetic held in which 
they diffuse. The relativistic electrons are thought to be 
accelerated in supernova remnants, which result from the 
massive stars. Given the many disparate steps, the com¬ 
pletely different physics, and the different time scales for 
the various processes involved, the tightness of the global 
FIR-RC correlation is surprising indeed. 

What is known about the relationship of molecular gas 
to the FIR and RC distributions? The FIR emission is well 
correlated with the CO emission on large scales in galaxies 
(e.g. Devereux & Young 1990; Young & Scoville 1991 and 
references therein), and the CO emission is also well cor- 
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related with the RC on global scales (Rickard, Turner & 
Palmer 1977; Israel & Rowan-Robinson 1984; Adler, Allen 
& Lo 1991; Murgia et al. 2002). At intermediate resolu¬ 
tions of ^1-3 kpc. Marsh & Helou (1995) found a small 
but systematic decrease in the FIR/RC ratio as a function 
of increasing radial distance, which they interpreted as a 
smearing of the RC by the propagation of the radiating 
electrons. By contrast, Adler et al. (1991) and Murgia et 
al. (2002) found that the ratio between the RC and the 
CO emission is constant, to within a factor of 3, both in¬ 
side the same galaxy and from galaxy to galaxy, down to 
kpc size scales. 

The arcminute-scale resolution of IRAS and ISO has 
limited extragalactic comparisons of the FIR-RC correla¬ 
tion to kpc sizes for all but a few Local Group sources. 
However, on size scales of a few hundred parsecs and 
smaller, detailed maps of gas within several kpc of the Sun 
suggest that the radio continuum emission associated with 
regions of star formation is, for the most part, thermal, 
and does not follow the FIR-RC correlation (Boulanger 
& Perault 1988; Haslam & Osbourne 1987; Wells 1997). 
Why there should be a good global correlation between 
the CO, FIR and the synchrotron-dominated RC is thus 
not at all expected from observations of the local interstel¬ 
lar medium. Furthermore, given the good correlations on 
large scales in galaxies, and the poor local correlation, one 
might expect that there may be some linear size scale(s) at 
which the CO-RC and FIR-RC correlations break down. 

In this paper, we examine several aspects of the CO- 
RC-FIR correlation problem. First, we quantify the FIR- 
RC, CO-RC and CO-FIR correlations on global scales for 
a sample of 24 galaxies with improved integrated single¬ 
dish CO measurements. Then, using &' imaging for a sub¬ 
sample of 9 sources, we show that the CO-RC relation is 
equally tight within galaxies as it is from galaxy to galaxy, 
down to linear scales of ~100 pc. We show that the most 
significant large-scale variations in the CO-RC relation ap¬ 
pear to be associated with spiral arms and that there is not 
a significant radial variation in the ratio. Finally, we pro¬ 
pose that the RC-FIR-CO correlations may be explained 
on large and small scales to be a result of regulation by 
the hydrostatic pressure in galaxies. 


2. Sample Selection and Data Reduction 

Galaxies were selected from the Berkeley-Illinois- 
Maryland Association Survey of Nearby Galaxies (BIMA 
SONG), an imaging survey of the CO (J=l-0) emission 
in 44 nearby spiral galaxies with a typical resolution of 
6" (Regan et al. 2001; Heifer et al. 2003). For 24 of the 
galaxies in the BIMA SONG sample, fully-sampled single¬ 
dish CO images were incorporated into the BIMA SONG 
images, so that all of the CO flux is imaged for these 
galaxies (i.e., there is no “missing flux” problem for these 
sources). The single-dish data for these 24 sources, taken 


at 55" resolution with the NRAO 12m telescope^, com¬ 
prise the most extensive collection of fully-sampled, two- 
dimensional single-dish CO images of external galaxies to 
date, and we use global CO flux measurements from these 
single-dish data. 

From the BIMA SONG data set, we selected a sub¬ 
sample of nine CO-bright galaxies for which high qual¬ 
ity VLA B-array at 1.4 GHz were available either from 
the NRAO VLA archive^ or from the Faint Images of 
the Radio Sky at Twenty-Centimeters (FIRST) survey 
(White et al. 1997). The VLA archival data were cali¬ 
brated and imaged with the NRAO package AIPS follow¬ 
ing standard procedures. To improve the sensitivity of the 
VLA radio continuum images to large scale structures, D 
array data were combined with the B array data. The D 
array data were taken from the VLA archive when avail¬ 
able, and combined in the visibility plane with the B array 
data; otherwise D array images from the NRAO VLA Sky 
Survey (NVSS, Condon et al. 1998) were added in the 
image plane using the linear combination method imple¬ 
mented in the MIRIAD task IMMERGE (Sault, Teuben, 
& Wright 1995). 

The subsample of 9 galaxies considered in this work are 
listed in Tabled along with some of their physical param¬ 
eters. All of them have an image of the CO integrated in¬ 
tensity from the NRAO 12m telescope. We therefore have 
CO and RC images at two different angular resolutions: 
a “low” resolution (~ 55") pair obtained from the NRAO 
12m and VLA D-array, and a “high” resolution (~ 6") 
pair. The RC and CO image characteristics for the low 
and high resolution data sets are listed in Table El and 01 
respectively, with beam sizes listed in angular and linear 
scales. 

We emphasize that these data are among the most 
comprehensive and uniformly processed CO and RC data 
from normal galaxies. Furthermore, the angular resolu¬ 
tion of the CO and RC images are well matched, not just 
by filtering the data, but by having comparable intrinsic 
coverage over the relevant spatial frequencies. Both the 
low-resolution and high-resolution data sets had intrinsic 
beam areas where the RC and CO images matched typi¬ 
cally to within 10%. Nonetheless, to account for even small 
differences in the beam areas, we convolved the images to 
common angular resolution in the quantitative analysis 
below. 

Finally, in order to study the RC, CO and FIR global 
correlations, the integrated 60 /rm fluxes of the galaxies 
in our sample have been taken from the IRAS and IRAS 
HIRES (^ 1') databases (see Tabled] for references). 

3. Data Analysis 

In order to compare the point-to-point RC and CO bright¬ 
nesses across entire galaxy disks, we overlaid regular grids 

^ The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated nnder cooperative 
agreement by Associated Universities, Inc. 
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Table 1. Galaxy Sample 


Galaxy 

a(J2000) 

{h m s) 

(5(J2000) 

(O , 

d 

(Mpc) 

(°) 

P.A. 

(°) 

RG3 Type, Nucleus 

IC 342 

03 46 49.7 

68 05 45 

3.9 

31 

37 

SAB(rs)cd, HII 

NGG 4258 

12 18 57.5 

47 18 14 

8.1 

65 

176 

SAB(s)bc, Syl.9 

NGG 4414 

12 26 27.2 

31 13 24 

19.1 

55 

159 

SA(rs)c?, T2: 

NGG 4736 

12 50 53.06 

41 07 14 

4.3 

35 

100 

(R)SA(r)ab,T2 

NGG 5005 

13 10 56.23 

37 03 33 

21.3 

61 

65 

SAB(rs)bc, LI.9 

NGG 5033 

13 13 27.53 

36 35 38 

18.7 

62 

170 

SA(s)c, Syl.5 

NGG 5055 

13 15 49.25 

42 01 49 

7.2 

55 

105 

SA(rs)bc, T2 

NGG 5194 

13 29 52.35 

47 11 54 

7.7 

15 

0 

SA(s)bc pec, Sy2 

NGG 6946 

20 34 52.33 

60 09 14 

5.5 

54 

65 

SAB(rs)cd,HII 


Coordinates, distances, inclinations, position angles, RC3 type and nuclear classifications from Heifer et al. (2003). 


Table 2. Low resolution CO and RC image characteristics. 


Galaxy 

VLA data 

RC Omai X Omin 
("X") 

a RC 

(mjy bm”'^) 

CO Ornaj ^ ^min 
("X") 

n ■ X D 
(kpc X kpc) 

b) 

'^CO 

(Jy bm“^ km s“^) 

IG 342 

D 

52x44 

0.15 

55x55 

1.0 X 1.0 

23.3 

NGG 4258 

D 

45x45 

0.50 

55x55 

2.1 X 2.1 

31.4 

NGG 4414 

C 

14x14 

0.20 

55x55 

5.0 X 5.0 

36.2 

NGG 4736 

D 

54x54 

0.13 

55x55 

1.1 X 1.1 

18.4 

NGG 5005 

D 

48x48 

0.20 

55x55 

5.6 X 5.6 

54.7 

NGG 5033 

D 

54x54 

0.20 

55x55 

4.9 X 4.9 

44.8 

NGG 5055 

D 

54x54 

0.20 

55x55 

1.9 X 1.9 

35.0 

NGG 5194 

D 

54x54 

0.30 

55x55 

2.0 X 2.0 

24.7 

NGG 6946 

D 

44x40 

0.10 

55x55 

1.5 X 1.5 

35.0 


a) The linear resolution refers to the CO beam. 

b) cco is the rms level in the undipped image of integrated intensity. 


of rectangular beam-sized boxes on both the RC and CO 
images. We then calculated the 1.4 GHz RC brightness 
(Iia) and the CO integrated intensity (/co) by averag¬ 
ing all pixel values within each box. For the 55" reso¬ 
lution data, we used CO images of integrated intensity 
made without clipping any of the channel images. For the 
6" resolution images, we used CO images as presented in 
Heifer et al. (2003), made by applying a smooth-and-mask 
technique that is very effective at showing low-level emis¬ 
sion (but that may introduce a bias against faint, com¬ 
pact emission; see van Gorkom & Ekers 1989; Heifer et al. 
2003). The smooth-and-mask technique does not bias the 
noise statistics of the final image, and we used the mask 
to calculate the pixel-by-pixel rms in the high-resolution 
images. For the radial profiles (below), the signal to noise 
ratio of the high- and low-resolution CO data is sufficient 
to use the undipped images of the integrated intensity, 
avoiding any bias. 

Fig. 1 shows the data analysis performed for the nine 
galaxies. For each galaxy, the upper and bottom left panels 
show the overlay of Iia (contours) on Ico (color) images 


respectively at low and high resolution. On top of the two 
panels, the intensity scale for the CO images is shown as 
a wedge. The RC contours start from a level of 5o- and 
scale by a factor of v^. 

The upper right panels plot Iia versus Ico at low 
(solid dots) and high (open dots) resolution. In the bot¬ 
tom right corner insets, the upper and lower panels show 
the box grids for the low and high resolution images, re¬ 
spectively. Only points above 2a have been plotted. Areas 
with obvious background continuum sources were omitted 
from the analysis. 

Finally, the bottom right panel shows the radial pro¬ 
files of IiA and Ico (on the top) and their ratio Qco/rc= 
-^Co(Jy sr“^km s“^)//i, 4 (mJy sr“^) (on the bottom). For 
this analysis, we used undipped images of integrated in¬ 
tensity in CO to avoid any bias against faint emission, 
and we calculated the mean intensities in concentric el¬ 
liptical annuli using the position angles and inclinations 
listed in Table Q We spaced the annuli to have widths 
along their major axis at the resolution of the correspond¬ 
ing images, or 6" for the high-resolution images and 55" 
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Table 3. High resolution CO and RC image characteristics. 


Galaxy 

VLA data 

RC Omaj X ^min 

{"x") 

O'RC 

(mjy bm“^) 

CO Omaj X Omin 

n y n 

(pcxpc) 

b) 

'^CO 

(Jy bm“^ km s“^) 

IC 342 

B-bD 

5.5X5.0 

0.05 

5.6X5.1 

110x97 

4.8 

NGG 4258 

AnB+D 

3.5X3.3 

0.05 

6.1x5.4 

240x210 

4.8 

NGG 4414 

B-bD 

5.4X5.4 

0.08 

6.4x5.0 

590x460 

1.9 

NGG 4736 

B-bD 

6.1X5.1 

0.05 

6.9X5.0 

140x100 

3.5 

NGG 5005 

B-bD 

5.4X5.0 

0.03 

6.2x6.0 

640x620 

2.9 

NGG 5033 

B-bD 

5.3X5.0 

0.07 

6.1x5.4 

550x490 

5.6 

NGG 5055 

B-bD 

5.5X5.0 

0.04 

5.8X5.5 

200x190 

3.6 

NGG 5194 

B-bG 

6.0x4.8 

0.02 

5.8X5.1 

214x188 

3.1 

NGG 6946 

B-bG 

5.6x4.9 

0.10 

6.0X4.9 

160x130 

3.1 


a) The linear resolution refers to the CO beam. 

b) crco is the rms level in the undipped image of integrated intensity. 


for the low-resolution images. We plot only those points 
that are measured at a level >la. For each galaxy, the val¬ 
ues of the mean and dispersion of Qco/RC= log((5co/Rc) 
are given on the bottom right plot of the radial variation 
in Oco/RC- 


4. Results 

4.1. Comparison of Global Fluxes 

For the 24 BIMA SONG galaxies that have associated 
single-dish CO images, Heifer et al. (2003) tabulated CO 
global fluxes that were directly measured from the single¬ 
dish images. For these sources, we can quantify the global 
CO-RC, FIR-CO, and RC-FIR correlations. The relevant 
fluxes are listed in TableQ] and the relations are plotted in 
Fig.ia Linear least squares fits to the relations yield a CO- 
RC correlation of log(Si, 4 )/log(Sco) oc 1.28±0.13, a CO¬ 
FIR correlation log(SF/fi)/log(Sco) cx 1.05±0.11, and a 
FIR-RC correlation of log(Si, 4 )/log(S_F/R) oc 1.05±0.14. 
(NGC 1068 is omitted from these fits, since it has a sig¬ 
nificant contribution to its RC flux from its AGN; how¬ 
ever, omitting this source has a small impact on the slopes 
derived.) The RC-FlR correlation measured for these 23 
sources is consistent with its best determination over a 
large (1809) number of sources (Yun, Reddy & Gondon 
2001), where the slope of the RG/FIR fit (also using 1.4 
GHz RC and 60 ^m FIR) is 0.99±0.01. 

Of the three correlations measured here, the slope of 
the CO-RC correlation is the steepest, and is the only 
one that appears to be inconsistent with a slope of unity. 
Furthermore, the CO-RC slope is consistent with that 
measured by Murgia et al. (2002), or 1.3 ± 0.1, on ar- 
cminute scales using independent data. However, given the 
uncertainties, it is fair to say that these data show that all 
three correlations are comparably tight, and that we can¬ 
not identify any of the relations as being a “fundamental” 
one. 


4.2. Polnt-to-Point and Radial Correlations of the RC 
and CO Emissions 

The 55" RC and CO images are shown in the upper-left 
panels of Figure 1. The RC emission, shown as contours in 
these figures, occasionally includes point sources which are 
almost certainly background radio sources. Ignoring these 
background sources, we can already see by eye that the 
general shape of the RC contours follows the underlying 
CO emission reasonably well. There are a few exceptions, 
like NGC 4258, for which the RC and CO emission seem 
to display quite different behaviours (see the Appendix for 
details). 

The 6" resolution images are shown in the lower-left 
panels of Figure 1. Here, the detailed agreement is even 
more striking, partly because a larger dynamic range of 
emission is represented in these images. We note that both 
the high- and low-resolution RC images have a larger dy¬ 
namic range than the corresponding CO images; this pre¬ 
vents a full comparison of the two in the faintest regions 
and it may give the impression that there are RC features 
that are unmatched by the CO. However, the opposite 
is not observed: every CO feature appears to have a RC 
counterpart. 

The quantitative point-to-point comparison of the 
/i ,4 — Ico correlation, shown for both the low and high 
resolution images in the upper-right panels of Figure 1, 
shows an excellent agreement between the two data sets 
for each galaxy. First, we note that for the point-to-point 
/i .4 — Jco correlation, the low-resolution data are aligned 
with the extrapolation of the high-resolution data to lower 
intensity values. Moreover, the /i ,4 and Ico radial profiles 
of the low and high resolution data also match very well 
(see bottom-right panels of Fig. 1). This is particularly ev¬ 
ident toward the central region of the galaxies where the 
two data sets overlap and is an effect that results from the 
incorporation of extended flux from the low-resolution im¬ 
ages into the high-resolution data. 
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Fig. 1. Radio continuum and CO comparison in BIMA SONG galaxies. (Upper left) RC contours overlaid on CO 
emission for 55" images. (Lower left) RC contours overlaid on CO emission for 6" images. Contours start at 5cr and 
scale by a factor of v^; the bright contours are chosen to match the dynamic range of the CO emission. (Upper right) 
IiA as a function of Iqo for the 55" and 6" data sets. The insets show the box grids used; the grid size corresponds to 
the resolution, so that each point is essentially independent. The dotted (dashed) lines represent the mean (dispersion) 
of the RC/CO ratio for all nine sources. The solid blue line is a weighted fit to the points shown which takes into 
account the errors in both coordinates. (Lower right) Radial profiles of /co and /i ,4 (top) and Qco/rc — Ico/h.A 
(bottom) for the 55" and 6" images. The mean and dispersion of Q'co/RC — fog(Qco/Rc) is also given. See §3 for 
further details on the presentation; notes on individual galaxies are given in Appendix A. (a) IC 342. 


The mean and dispersion in qco/KC are given in the 
bottom right panels of Figure 1. The dispersion within the 
individual sources ranges from 0.09 to 0.35; these corre¬ 
spond to dispersions in the CO/RC flux ratios that are a 
factor of 1.2-2.2. These are remarkably tight correlations 


given that the dynamic range extends over three orders of 
magnitude for some sources. The average gco/RC among 
all nine sources is 1.03 ± 0.06 (uncertainty in the mean) ± 
0.17 (standard deviation); that is, the CO/RC ratios from 
all sources are constant to within a dispersion of a fac- 
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Fig. 1. (b) NGC 4258. See Figure la. 

tor of 1.5. These results emphasize both the small scatter 
in qco/RC within each galaxy as well as among different 
sources. The CO and RC emission are on average linearly 
related for these sources, though we note that larger sam¬ 
ples show a steeper RC/CO slope (§ 4.1). 

Despite small-scale variations, there is no systematic 
trend of qco/RC with radius (lower-right panels of Figure 
1). In FigureEl we show the correlations of all of the galax¬ 
ies on a common plot. The individual points represent an¬ 
nular averages and may generally be read as radial profiles. 
This figure serves to reinforce the commonality of the CO- 
RC relation among the galaxies. We note that this figure 
shows distance-independent intensity measurements, not 


atopw-'OO? rrTJ«*0 33 d«x 



(Udiun (arcsec) 

0 50 100 ISO 



luminosity measurements, and that the intensities extend 
over 3 orders of magnitude. 

Figure 121 shows the correlation in the global (7co/RC 
over entire galaxies. Figure O shows the correlation at in¬ 
termediate scales, when averaged over annular bins. We 
can examine the scatter in qco/RC on the smallest scales 
by placing the individual pixels from all sources on a com¬ 
mon plot, which we present in Figure 0] Over all pixels, 
the average correlation is <9co/RC> = 1-1 i 0.28. 

Some of the galaxies show slightly depressed CO emis¬ 
sion relative to the RC at their centers, which may be 
caused by enhanced RC emission from the presence of an 
ACN. But with the exception of IC 342, the deepest de- 
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Fig. 1. (c) NGC 4414. See Figure la. 

pression in the emission ratios is only a factor of 1.1-3.3 
away from the mean for each source. This close agree¬ 
ment in the galaxy centers is somewhat surprising, given 
the likely variation in the CO-to-H 2 conversion in these 
regions. 

To summarize, we find a very good point-to-point cor¬ 
relation between the CO and RC emissions at an angluar 
resolution of 6" in a sample of 9 nearby spiral galaxies. We 
are aware that the small size of the sample analyzed in this 
work could weaken the statistical relevance of this result. 
However, for a sample of 22 galaxies that includes an ad¬ 
ditional 13 sources, Paladino et al. (in preparation) find a 
similarly low scatter in the overall RC-CO correlation. 


8topw"'0 67 rrTJ8*0Jl drx 



(UdiUA (arcsftc) 

0 50 100 



HadluB ikpc) 


4.3. Variations in qqqj-^q With Spiral Structure 

To investigate trends as a function of location within the 
galaxies, we also made images of the ratios qco/RC for 
each source (FigureEJ. All images are shown on a common 
color scale, so that color variations from source to source 
indicate real variations in qqq . 

There are clearly some differences in the distribution 
of Qco/RC among the different sources. For NGC 4414, 
NGC 5005, and NGC 5033, there is very little variation in 
Qco/rCj as evidenced by their low dispersions (Fig. 1). It 
is worth noting that both NGC 4414 and NGC 5005 are 
flocculent galaxies (Elmegreen & Elmegreen 1987), which 
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Fig. 1. (d) NGC 4736. See Figure la. 


lack prominent spiral structure. In NGC 4736 and NGC 
5194 (M 51), systematic large-scale variations are appar¬ 
ent that follow the galaxies’ spiral structure. This is seen 
in NGC 6946 as well, with enhanced 9co/RC along the 
central north-south elongated structure (variously inter¬ 
preted as a bar or as inner spiral arms; see Regan & Vogel 
1995) and in the somewhat wispy arms to the southwest. 
In NGC 5055, there is a large-scale ringlike enhancement 
in Qco/Rc; remarkably, the enhancement in this floccu- 
lent prototype appears to follow nearly exactly two inner 
spiral arms revealed in K' imaging (Thornley 1996). The 
most significant, large-scale structures in 9co/RC therefore 
appear to be organized along spiral arms. 


The two galaxies with the highest average qqq /rc 
among the highest dispersion in gco/RC IC 342 and 
NGC 4258. The variations in qco/RC in these galaxies do 
not appear to be closely correlated with spiral structure; 
however, these galaxies are unusual for other reasons, dis¬ 
cussed further in Appendix A. 

The association of higher CO/RC ratios on 1-2 kpc 
scales in spiral arms was noted by Adler et al. (1991), 
who also measured a larger scatter in the CO/RC rela¬ 
tionship in these regions. Given the improved resolution 
of our images, we performed a further analysis of the RC- 
CO correlation behavior in the spiral arms of M51 with 
the intent of determining whether this systematic varia- 
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Fig. 1. (e) NGC 5005. See Figure la. 


tion in 9co/RC might correspond to a breakdown in the 
correlation. We emphasize that we are sampling the RC 
and CO intensities at a spatial resolution of about 200 pc, 
a value which is significantly below the electron diffusion 
scale length. The result of this investigation is shown in 
Fig.El M51 still shows a tight linear RC-CO correlation, 
with a scatter of only 0.2 dex, within the arms. Based 
on this result, we cannot conclude that correlation breaks 
down at the linear scale of an arm width. If we consider 
the distribution of the CO/RC ratios in the arms, we get 
9co /RC = 1.2 ± 0.2, which is only slightly above the aver¬ 
age of 1.1 ± 0.28 that we get for the entire galaxy sample. 
In the interarm regions, we measure qco/RC =1-0 ± 0.25. 


On top of this systematic effect, which appears to be asso¬ 
ciated with spiral structure, we observe larger but isolated 
deviations in 9co/RC which contribute to the correlation 
scatter. However, despite these localized variations there 
is still a remarkably low dispersion in the high-resolution 
CO-RC correlation. To summarize, we tend to measure 
an increase in CO/RC in spiral arms as seen earlier in 
Adler et al. (1991); unlike Adler et ah, we do not see an 
increase in scatter within the spiral arms: the correlation 
dispersion is just as low. 















10 


M. Murgia et al.: The Molecular Connection to the FIR-Radio Continuum Correlation in Galaxies 



200 400 «0D 100 

JV.fiEAMKMfi 



200 lOD 0 *100 >200 

ReUT1\ERA IARCS6C1 



10 “ 


^ ,0.0 

n 

u 

b 

> 10- 

i 

J ,o» 


10 * ,o»* io>* 10 “ 

(Jy/olernd km/*) 



(Udtu* (nrcftcc) 

0 too 200 300 



Fig. 1. (/) NGC 5033. See Figure la. 


4.4. The CO, RC and FIR Length Scales 

We have modeled the CO and RC exponential scale 
lengths from 22 BIMA SONC galaxies with CO data taken 
from the 55" 12 m data and the low-resolution RC data 
from either the present study or from NVSS imaging. The 
results are listed in Table 0] and shown in Figure [Tj We 
note that the formal errors in the scale lengths almost cer¬ 
tainly underestimate the true errors, for the following rea¬ 
sons: (1) the images often cover only 1-2 e-folding lengths 
in the galaxies; (2) given the small coverage, the radial 
profiles may or may not be fairly represented by exponen¬ 
tial functions. At higher resolution, the CO radial profiles 


have much more detailed structure and are in general not 
well characterized by exponentials (Regan et al. 2001); (3) 
there may be systematic errors in the CO data that are 
not well represented by the formal errors (Heifer et al. 
2003); (4) we have not excluded background radio con¬ 
tinuum sources; this error might tend to produce longer 
RC scale lengths than are present intrinsically. We note 
that some of the CO scale lengths derived here differ from 
those published in Regan et al. (2001). The scale lengths 
presented here are to be preferred (Regan 2004, private 
communication). 

Despite the qualifications (which may also apply to 
previously determined FIR and RC scale lengths), the 
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Fig. 1. {g) NGC 5055. See Figure la. 


data in Table 0 suggest that on average, the CO scale 
lengths are comparable to the RC scale lengths, with 
< Ico/Irc > = 0-98 ± 0.41. This is true even though 
the dispersion in the average is large; that is, some galax¬ 
ies show significantly longer or shorter CO disks than RC 
disks. We can further compare the FIR scale lengths mod¬ 
eled by Marsh & Helou (1995) with the BIMA SONG CO 
disk scale lengths for an overlap sample of 5 galaxies: NGC 
2903, NGC 4303, NGC 5055, NGC 5194, and NGC 6946. 
(We exclude NGC 628, which lacks RC in the central re¬ 
gion where the CO is present; the fit to the ringlike RC 
distribution yields a large scale length which is poorly de¬ 
termined. We also exclude NGC 7331, with ringlike emis¬ 


sion in CO and RC that also precludes a good determina¬ 
tion of its scale length.) After scaling the Marsh & Helou 
(1995) leofim to a common distance scale with the BIMA 
SONG study, we find that for this subsample, l^o^m/lco 
= 0.69 ± 0.07, leofim/lRC = 0.65 ± 0.14, and Ico/Irc = 
0.94 ± 0.20. These results suggest that the FIR tends to 
have a smaller scale length than the RG or GO and that 
the RG and the GO may in some sense be more directly 
related to each other than either is to the FIR. New deter¬ 
minations of the FIR scale lengths using Spitzer imaging 
should help to clarify this issue. 
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Fig. 1. (h) NGC 5194. See Figure la. 


5. Discussion 

5.1. The Radial Correlations 

The most significant deviation from the close correlation 
between the FIR and RC emission is a small, generally 
monotonic decrease in FIR/RC as a function of radius in 
many extragalactic sources (Beck & Golla 1988; Rice et al. 
1990; Bicay & Helou 1990). In a seminal study, Marsh & 
Helou (1995) studied fully 2-dimensional maps of the FIR 
and RG emissions from 25 nearby galaxies, with improved, 
^ I' resolution over previous studies. Their results show 
decreasing ratios of FIR/RG in nearly all sources, with 
disk values at ~4' radius (10-20 kpc) that are lower by 


up to a factor of 10 than the central ratios. In the Marsh 
& Helou study, there are 3 galaxies in common with the 
present study: NGC 5055, NGC 5194, and NGC 6946. 
These galaxies show monotonic decreases in FIR/RC of a 
factor of 3.7, 4.1, and 2.8 over 4-5' in radius. By contrast, 
the galaxies in the present study, including these three 
galaxies, do not show any apparent trend of CO/RC with 
radius (see bottom right panels of Figure 1). The CO-RC 
correlation therefore seems to be different from the FIR- 
RC correlation: the CO/RC shows no systematic trend 
with radius, and the CO and RC scalelengths seem to be 
on average equal, with both longer than the FIR scale- 
lengths (§ 4.4). 
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Fig. 1. (i) NGC 6946. See Figure la. 


In the CO-RC correlation, we have seen that the 
strongest organized structures in qqq are found along 
spiral arms, with enhanced qco/RC along the spiral arms 
and lower values in the interarm regions. Recent Spitzer 
FIR imaging of M 81 (Gordon et al. 2004) shows very 
similar features: first, that the variations in (7 fir/rc are 
organized on size scales far greater than the resolution of 
the observations; and second, that these large-scale vari¬ 
ations in gpiR/RC strongly follow the spiral structure in 
M 81. In fact, though Gordon et al. don’t discuss this is¬ 
sue, an examination of the (?fir/rc map of M 81 appears 
to show no systematic radial trend in (7 fir/rc- 


How can the M 81 FIR-RC result and the CO-RC re¬ 
sults presented in this study, namely that the strongest 
variations in the ratios are seen in azimuth across spiral 
arms, not as a radial variation, be justified with IRAS- and 
ISO-era observations that showed decreasing FIR-RC as a 
function of radius? If the strongest intrinsic variations are 
indeed across spiral arms, we note that the effect of these 
variations will generally be largest at small radius. That is, 
when averaged in azimuthal rings, the fraction of the total 
area in an annulus taken up by the spiral arms is much 
larger at small radius than it is at large radius. It may be 
that the measured radial falloff in FIR-RC is primarily a 
consequence of the low resolution of the FIR observations; 
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Fig. 2. Global CO-RC, FIR-CO, and RC-FIR fluxes for 24 BIMA SONG galaxies, shown with least squares fits (from 
which NGG 1068 is omitted). The data are listed in Table @| Grosses are galaxies included in the high-resolution 
analysis of this paper. The dashed-line represents a least squares fit restricted to the data of the 9 galaxies analyzed 
in this work and yields a GO-RG correlation of log(5'R,c) oc 1.06 ± 0.171og(S'co)- 



Fig. 4. Histogram of all measured qco/RC pixels. The 
mean and dispersion are also given. 


it is essentially a measurement of the arm-interarm con¬ 
trast of the galaxy. We predict that high-resolution Spitzer 
observations from the SINGS team (Kennicutt et al. 2003) 
and others will show that M 81 is not anomalous in this 
sense: the dominant variation will also be across spiral 
arms, with any radial trend being a secondary effect. In 
any case, Spitzer imaging will certainly further illuminate 
the issue of the relative strengths of the GO-RG, FIR-RG 
and the GO-FIR correlations. 


5.2. A Brief Summary of Models of the FIR-RC and 
CO-RC Correlations 

The global FIR-RG correlation is valid over five orders 
of magnitude in luminosity (Yun et al. 2001). We have 
shown in this paper that the RG, FIR and GO emission 
are also globally correlated, and that the GO-RG corre¬ 
lation is equally good on scales down to several hundred 
parsecs. However, the physical bases for understanding the 
FIR-RC and CO-RC correlations are not well understood. 
The problem is that the emission mechanisms are so dif¬ 
ferent that it is difhcult to understand why such a tight 
relationship should hold between them. The FIR emis¬ 
sion is thought to be from dust heated by UV radiation 
from hot young stars and is therefore a reasonably good 
tracer of the star formation rate in a galaxy (e.g. Devereux 
& Young 1991). Because the extinction in the UV is so 
high, essentially all of the UV photons are absorbed by 
dust and are reprocessed into FIR emission. Few Lyman 
continuum photons are thought to escape from a galaxy 
(Leitherer et al. 1995; Heckman et al. 2001; Deharveng et 
al. 2001). Over scales of a few hundred parsecs or more, the 
RC emission from galaxies is typically at least 90% non- 
thermal (Condon 1992); for the galaxies in this paper, the 
thermal fractions at 1 GHz are typically 5-10% (Niklas, 
Klein & Wielebinski 1997). This synchrotron emission is 
due to cosmic-ray electrons spiraling around the galactic 
magnetic field. Why should the emission from relativistic 
cosmic-ray electrons be so tightly correlated with the FIR 
from dust? The CO emission from giant, turbulent, cloud 
complexes traces the bulk of molecular gas mass, even 
though it is optically thick. The good correlation between 
the CO and the FIR has been known for some time and 
is generally attributed to the intimate connection between 
massive star formation and molecular clouds (Devereux & 
Young 1991). That is, because CO has been found to be a 
reliable tracer of the mass of GMCs, and because the rate 
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slope=1.08 rms = 0.27 dex 



Fig. 3. CO-RC correlation in sample galaxies. Each point 
line is a weighted fit to the points shown which takes into 

1.08. 

of massive star formation is related to the mass of molecu¬ 
lar gas through an approximately constant star formation 
efficiency averaged over time (Rownd & Young 1999), the 
FIR, through the UV production and star formation rates, 
should be coupled to the CO emission. In detail, however, 
there is scatter of at least an order of magnitude in the 
CO-FIR relation (Young & Scoville 1991), and it has been 
difficult to relate the FIR reliably to a star formation rate 
for different types of galaxies (e.g. Kennicutt 1998; Bell 
2003). 

In the following, we briefly review some of the mod¬ 
els previously posed to explain the FIR-RC and CO-RC 
correlations. 

5.2.1. The FIR-RC Correlation: Conventional Picture 
and Calorimetric Models 

In the conventional model, the FIR and the RC are cor¬ 
related because each is a tracer of massive star formation: 
the RC through the supernovae that are ultimately pro¬ 
duced from the recently formed massive stars and the FIR 
from the dust-reprocessed UV radiation from the stars 
themselves. But this explanation is wanting on several ac¬ 
counts. First, why should the density of cosmic-ray elec¬ 


represents the average value in an annulus. The solid blue 
account the errors in both coordinates and has a slope of 


trons and the strength of the interstellar magnetic field be 
closely coupled to the rate of star formation? The produc¬ 
tion of the RC emission requires a series of steps: massive 
star formation —> supernova production —> cosmic-ray ac¬ 
celeration —!■ interaction with the magnetic field, as does 
the FIR emission: massive star formation ^ UV photon 
production ^ absorption and reradiation by dust. To get 
the FIR-RC correlation to hold to a factor of two, as ob¬ 
served, each step must have a smaller range of variation 
than the observed correlation, which seems implausible. 
Furthermore, since the FIR-RC correlation is observed 
over a large range of physical scales, it is unclear why 
it should be so good on galactic scales, since star for¬ 
mation is a local process. In addition, star formation is 
seen to correlate well with thermal RC emission in HII 
regions (Boulanger & Perault 1988; Haslam & Osbourne 
1987; Wells 1997), rather than the non-thermal emission 
on scales less than 100 pc, and there is not a particularly 
good correlation of supernova remnants with either sites 
of star formation (Wootten 1978) or with mid-IR or FIR 
emission (Cohen & Green 2001; Whiteoak & Green 1996). 

Another problem with the conventional model is that 
the CO and RC both have a somewhat poorer correla¬ 
tion with other SFR indicators than they do with each 
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Fig. 5. Spatially resolved CO-RC correlation within NGC 5194 spiral arms. The RC and CO intensities are measured 
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1.04 and with a scatter which is less than a factor of 2 . The RC/CO ratio in the arms is 9co/RC= 1-2 ± 0.20, while 
in the interarm regions (Zco/RC= 1-0 ^ 0.25. 


other; indeed Bell (2003) raises issues about the linearity 
of the FIR-RC correlation by arguing that both indicators 
underestimate the star formation rates at low luminosi¬ 
ties. Bell goes on to reach the unlikely conclusion that two 
unrelated effects have conspired to generate the observed 
correlation, without questioning the fundamental assump¬ 
tion that it is the star formation rate that provides the key 
link to the FIR-RC correlation. Leroy et al. (2004) show 
that the CO-RC correlation holds even for low-luminosity 
dwarf galaxies, suggesting that even further ad-hoc tuning 
would be necessary in order to understand the correlations 
among the CO and the FIR and RC within a conven¬ 
tional SFR framework. Notwithstanding Bell’s arguments 
regarding the suitability of the FIR and RC as star for¬ 
mation tracers, the correlations among the FIR-RC and 
CO-RC remain tight even if the FIR and RC underesti¬ 
mate the star formation, and we must look to an alternate 
explanation to understand the observed correlations. 

Volk (1989) and others since have presented a calori¬ 
metric model in which the ratio of RC and FIR emissions is 
held constant by assuming that both the relativistic elec¬ 
trons as well as the UV photons are proportional to the 
SN rate. In this scenario, all systems are optically thick, 
so that all synchrotron or inverse Compton losses must 
occur in the host galaxy, and all the FUV radiation from 
massive stars is absorbed by dust grains within the galaxy. 
With both of these processes driven by the massive stars, 
the galaxy acts as a calorimeter. 


Niklas & Beck (1997) show that only ~ 30% of their 
74-galaxy sample have radio spectral indices that sup¬ 
port the optically thick model; however, Lisenfeld & Volk 
(2000) argue that cosmic ray electron convection, rather 
than diffusion, can flatten the spectra within the calorime¬ 
ter model. Finally, the calorimeter model addresses the 
global energetics within galaxies, and it does not address 
the observed FIR-RC correlation that holds at least on 
kiloparsec size scales within galaxies, or the CO-RC cor¬ 
relation that holds down to size scales of ^ 100 pc. This 
model thus seems to be unsatisfactory as a solution that 
holds generally for different kinds of galaxies at various 
resolutions. 

In the phenomenological model by Helou & Bicay 
(1993), the FIR-RC correlation is driven by a tight cou¬ 
pling between the dust-heating photon luminosity and the 
production rate of cosmic ray electrons. The photon lumi¬ 
nosity is regulated by the effective optical depth to the 
photons, and the CR luminosity is regulated by the effi¬ 
ciency with which the cosmic rays emit synchrotron radi¬ 
ation before they escape. In the optically thick case, the 
model is similar to calorimetric models, which Helou & 
Bicay postulate applies to infrared-luminous galaxies. In 
the optically thin case, which Helou & Bicay apply to 
normal galaxies, the photon luminosity and the CR lumi¬ 
nosity are predicted to be approximately equal because 
of a coupling of the magnetic field with the gas density. 
Niklas & Beck (1997) however point out that in this model, 
the transition between galaxies with low and high radia- 
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Table 4. Global Fluxes and Scale Lengths in BIMA SONG Sources 


Galaxy 

S1.4 

(mJy) 

Seo 

(Jy) 

Sco 

(Jy km s“^) 

Reference 

‘‘CO 

(kpc) 

I 

‘‘RC 

(kpc) 

Ico/Irc 

IC 342 

2475 

255.96 

20000 

1,2,3 

3.5 ± 0.2 

6.21 ± 0.06 

0.56 

NGG 0628 

180 

20.86 

1514 

2 

3.2 ± 0.1 

15.3 ± 0.1 

0.21 

NGG 1068 

4991 

181.95 

4102 


1.6 ± 0.4 

0.7 ± 0.2 

2.16 

NGG 2903 

407 

60.03 

3254 


2.0 ± 0.2 

1.7 ± 0.2 

1.21 

NGG 3351 

48 

19.92 

1513 


- 

- 

- 

NGG 3521 

357 

47.02 

4800 


2.1 ± 0.2 

4.1 ± 0.1 

0.52 

NGG 3627 

458 

67.8 

4259 


5.1 ± 0.1 

3.5 ± 0.1 

1.45 

NGG 3938 

62 

9.24 

923 


4.0 ± 0.2 

6.3 ± 0.1 

0.64 

NGG 4258 

792 

21.60 

2686 

1,2 

2.1 ± 0.2 

3.1 ± 0.3 

0.70 

NGG 4303 

416 

37.53 

2427 


3.2 ± 0.1 

3.3 ± 0.1 

0.96 

NGC 4321 

340 

25.86 

2972 


4.8 ± 0.1 

5.0 ± 0.1 

0.96 

NGC 4414 

231 

30.11 

2453 


3.4 ± 0.2 

3.00 ± 0.03 

1.13 

NGC 4569 

83 

10.08 

1096 


2.6 ± 0.1 

2.5 ± 0.2 

1.04 

NGC 4736 

254 

69.2 

2641 


0.9 ± 0.4 

0.7 ± 0.2 

1.22 

NGC 4826 

103 

35.45 

1845 


0.6 ± 1.1 

0.4 ± 0.7 

1.46 

NGC 5005 

176 

22.3 

1278 


2.9 ± 0.3 

3.2 ± 0.1 

0.90 

NGC 5033 

178 

16.45 

2469 


4.5 ± 0.1 

3.2 ± 0.1 

1.39 

NGC 5055 

390 

40.02 

3812 


2.4 ± 0.1 

2.3 ± 0.1 

1.03 

NGC 5194 

1490 

98.8 

10097 


2.9 ± 0.1 

3.4 ± 0.1 

0.84 

NGC 5247 

53 

8.07 

1102 

4,5 

7.1 ± 0.2 

7.3 ± 0.1 

0.97 

NGC 5248 

140 

20.71 

1829 


3.0 ± 0.3 

4.2 ± 0.1 

0.72 

NGC 5457 

750 

88.04 

3479 


- 

- 

- 

NGC 6946 

1395 

136.69 

9273 

2,4 

2.5 ± 0.2 

3.83 ± 0.04 

0.65 

NGC 7331 

373 

35.29 

2762 

2,4 

3.2 ± 0.1 

3.33 ± 0.03 

0.96 

Average 







0.98 ± 0.41 


а) Flux references: (1) Si ,4 is 1.4 GHz RG from White & Becker 1992; (2) Seo is 60 fim emission from Rice et al. 1988; (3) Sco 
from Crosthwaite et al. 2001; (4) Si ,4 from Condon 1987. (5) SeojSioo from IRAS Faint Source Catalog; If not otherwise 
indicated, Si .4 is from Condon et al. 1990, Seo is from the IRAS Bright Galaxy Sample, Soifer et al. 1989, and Sco is from 
Heifer et al. 2003. 

б ) GO scale lengths modeled from the 55” 12 m data. 

c) RG scale lengths modeled from the low-resolution RC data from either the present study or from NVSS imaging. 


tion fields should be accompanied by an increase in the 
nonthermal spectral index, which they do not observe. 
Furthermore, Niklas & Beck show that their measured 
spectral indices do not correlate well with the equipar- 
tition magnetic field strength, in contradiction with the 
Helou & Bicay (1993) picture. 

5.2.2. The CO-RC correlation: Cosmic Ray Heating 
and Secondary Electrons 

Suchkov, Allen & Heckman (1993) present a model for the 
GO-RG correlation where the nuclear component of cos¬ 
mic rays is responsible for heating the molecular gas while 
the electron component account for the synchrotron radi¬ 
ation. In this model, the CO need not be directly tracing 
star formation; rather, its emission depends on the cosmic 
rays that are products of high-mass star formation. This 
model is particularly attractive since it naturally provides 
a direct physical link between the CO line integrated in¬ 
tensity and the RC synchrotron emission. Bradford et al. 


(2003) show that the bulk of the molecular gas in the cen¬ 
tral 180 pc of NGC 253 is highly excited; these authors 
infer a temperature of T ~ 120 K for the molecular gas 
and conclude that the best mechanism for heating gas is 
cosmic rays that are produced towards the starburst nu¬ 
cleus of this galaxy. However, we show that the observed 
CO-RC correlation spans more than three orders of mag¬ 
nitude in both quantities. This poses a serious problem for 
the CR heating model in that the molecular gas temper¬ 
ature should be raised by at least a factor of 1000. Thus, 
if it cannot be excluded that the CR heating is acting 
in some circumstances, it is difficult to accept it as the 
general explanation for the CO-RC correlation. 

An interesting aspect of this picture, although not of¬ 
ten addressed in literature, is the issue of the secondary 
electrons, which are produced by CR primaries interact¬ 
ing with the molecular gas (e.g. Marscher & Brown 1978). 
If the production of the secondary electrons is significant, 
then the CR nuclei could not only raise the molecular gas 
temperature but also enhance the synchrotron radiation. 
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As outlined in § 1, the RC emission from isolated giant 
molecular clouds is mostly thermal. Thus, the secondary 
electrons able to survive collisional losses will probably 
escape molecular clouds in a very short time-scale and 
diffuse in the disk where they became indistinguishable 
from the primaries. One can devise a scenario in which 
the CO-RC correlation is governed by the volume density 
of molecular clouds. The CO integrated intensity traces 
the number of clouds within the volume of the galaxy disk 
intercepted by the beam. On other hand, the secondary 
electron injection will also be proportional to the volume 
density of clouds, leading to a linear CO-RC correlation. 
A precise calculation of the secondary electron contribu¬ 
tion implies special choices of rather uncertain parameters, 
such as the efficiency of energy transfer from supernovae 
to CRs and the ratio of primary electrons to nuclei in CRs 
(see e.g. Paglione et al. 1996) or the effective amount of 
diffuse molecular gas in the interstellar medium. While 
a detailed discussion of this issue is beyond the scope of 
this paper, such an intertwined relationship may help to 
explain the closeness of the CO and RC emission. 

5.2.3. Cosmic-Ray Driven Interstellar Chemistry 

Bettens et al. (1993) proposed that chemical evolution in 
molecular clouds, driven by the cosmic ray flux, is finked to 
the FIR flux through ion destruction in the clouds. While 
taking into consideration all three of the RC, FIR, and 
CO emissions, this theory neglects the strong influence 
of the interstellar magnetic field on the synchrotron ra¬ 
diation. The theory also makes rather specific predictions 
about how the relative abundances of molecules like HCN, 
H 2 CO, and CH 3 OH are tied to the nonthermal RC. Based 
on results from this paper and from how the HCN/CO ra¬ 
tios vary within the Milky Way and other galaxies (e.g. 
Heifer & Blitz 1997a;6), it seems unlikely that these pre¬ 
dictions will hold in detail. 

5.3. Hydrostatic Pressure as a Regulating Mechanism 
for the CO-RC Correlation 

The difficulties in relating the seemingly independent 
emission mechanisms, and doing so over the range of size 
scales observed, suggest that there might be a single phys¬ 
ical parameter that acts to correlate the observed quanti¬ 
ties, rather than some process that relates them directly. 
We suggest that the parameter is the midplane hydro¬ 
static pressure P, modified by other large scale pressure 
variations such as those produced by the spiral arms. 
Recently, Wong (2000), Wong & Blitz (2002) and Blitz 
& Rosolowsky (2004) have suggested that the azimuthally 
averaged fraction of gas in the molecular phase, fmoi, or 
equivalently, the H 2 to HI ratio, is determined by the ra¬ 
dial variation in interstellar gas pressure, which can be 
approximated as 

P = 0.84(GI]*)°%^£^- (1) 


Here Eg is the total surface density of gas, E* is the surface 
density of stars, Vg is the velocity dispersion of the gas, 
and /i* is the stellar scale height, which can be written in 
terms of the stellar velocity dispersion u* and midplane 
volume density p* as /i* = (u*^/ 47 rGp*)°-®. 

This expression for P assumes that the midplane vol¬ 
ume density of stars greatly exceeds the midplane volume 
density of gas, a condition that is satisfied throughout 
most disk galaxies except perhaps in spiral arms. One par¬ 
ticular utility of this formulation is that both Eg and E* 
are directly measurable, while Vg and /i* are observed to 
be fairly constant within galaxies {vg: Shostak & van der 
Kruit 1984; Dickey et al. 1990; /i*: van der Kruit & Searle 
1981a; 1981b; Fry et al. 1999). Furthermore, Vg shows lit¬ 
tle variation from galaxy to galaxy about a value of ~8 
km s“^, and /i* has a dispersion of only 50% from the 
mean (Kregel, van der Kruit & de Grijs 2002). In any 
event the dependence of the pressure on /i* is weak due to 
the exponent of 0.5. 

Let us first consider the origin of the CO-RC cor¬ 
relation. We begin by making the usual assumption 
that the CO surface brightness scales linearly with Eh 2 - 
Observations show that EHa/Eni oc _p0-8-0-9 (Wong & 
Blitz 2002; Rosolowsky & Blitz in preparation). Because 
the surface density of atomic hydrogen Ehi is roughly con¬ 
stant with radius (e.g. Wong & Blitz 2002), we find that 

Ico cx ® . ( 2 ) 

The shallow decrease of Ehi with radius in some galaxies 
would increase the exponent slightly. 

The RC emission, on the other hand, is determined 
largely by the synchrotron emissivity, which can be writ¬ 
ten as 

Jrc oc , (3) 

where B is the magnetic field strength, Nq is the num¬ 
ber density of cosmic-ray electrons with energies between 
energy e and e -I- de, and 7 is the power law index of the 
spectrum of cosmic-ray energies such that N(e) = NQe~'^ 
(Condon 1992). The power law index of the observed syn¬ 
chrotron spectrum suggests that typically 7 = 2.6 at GeV 
energies; thus we expect Jrc oc 

To express B in terms of P, we assume energy equipar- 
tition between the turbulent energy which dominates the 
hydrostatic pressure P and the magnetic energy density 
B^jSn, so that B cx .P° ®. This assumption has been dis¬ 
cussed previously in the literature with regard to the FIR- 
RC correlation (e.g. Helou & Bicay 1993; Groves et al. 
2003). We therefore expect Jrc oc NqP^-^. 

To express A^o in terms of P, we consider the radial 
distributions of both quantities as observed in the Milky 
Way. The radial distribution of Nq, although uncertain, 
is based on gamma-ray data from the Compton Gamma 
Ray Observatory, and yields a power-law fit: Nq oc 
for R ^ 3 kpc (Strong & Mattox 1996). The radial distri¬ 
bution of P is shown in Figure |H1 using the data of Dame 
(1993) for Eg and an exponential distribution of stars with 
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Fig. 6. Images of gco/RC as defined in §3. The scale is given as a color wedge at the top of the figure. 


a radial scale length of 3 kpc (Spergel, Malhotra & Blitz 
1996) tied to a value of S* = 35 Mq pc“^ at the solar 
circle. Beyond 4 kpc, P is seen to have an approximately 
power law form with P cx Combining these two 

relations suggests that Nq oc P^ "^^ , We therefore expect 

Ikc cx (4) 

Combining Equations 2 and 4 yields 

/rc cx Vo", (5) 

which is consistent within the uncertainties with the global 
results given in § 4.1 and with the intermediate-resolution 
results of Murgia et al. (2002) for spiral galaxies and of 
Leroy et al. (2004) for dwarf galaxies. 


As a check on the above derivation of iVo, we note 
that if we combine the above result Irq oc NqP^-^ with 
Equation (2), we derive Irc oc NqIqq. Eor the large- 
scale observations presented here, in Murgia et al. (2002) 
and in Leroy et al. (2004), we have /rc oc Iqo, which 
suggests a scaling of the CR electron density with pressure 
(A^o) cx (P) (quantities in brackets represent values 
averaged over the entire galaxy). For the sub-sample of the 
9 bright galaxies studied at high resolution in this work, we 
instead observe /rc cx Iqq- The sub-sample observations 
suggest that A^o rnay be nearly independent from P for 
these galaxies; in any case, the intrinsic variation from 
source to source suggests that A^g could vary somewhat 
among galaxies. However, given the small scatter of the 











20 


M. Murgia et al.: The Molecular Connection to the FIR-Radio Continuum Correlation in Galaxies 



Fig. 7. Histogram of the CO to RC scale lengths ratio. The 
average ratio considering all galaxies results 0.98 ± 0.41. 

CO-RC correlation, it appears that Nq varies only slowly 
with P. 

Next we consider how the FIR emission is related to 
both the CO and RC. The FIR emission has tradition¬ 
ally been linked to the CO emission through a Schmidt- 
type star formation law, with the assumption that the 
FIR is directly proportional to the star formation rate per 
unit area. This formulation works reasonably well for lu¬ 
minous FIR galaxies (Kennicutt 1998), but is problematic 
for normal and low-luminosity sources (Kennicutt 1998; 
Bell 2003). 

Recently Dopita et al. (2004) proposed that the FIR 
emission itself is also directly related to P: the ISM pres¬ 
sure sets the stall radius of the Hii region, so that with 
higher P, the H ii region stalls at smaller radius and the 
mean dust temperature in the atomic and molecular shell 
around the H ii region should be higher. From stellar spec¬ 
tral synthesis modeling, Dopita et al. show that the 60 
/im flux is positively correlated with the ISM pressure, 
though in their “toy” model the increase is not enough 
to account for the corresponding increase in the RC emis¬ 
sion. However, taken together with the above model of the 
CO-RC correlation, these models may have the potential 
to link all three of the CO-RC-FIR emissions with the in¬ 
terstellar pressure, without any explicit dependence on a 
Schmidt-type law or other star formation scenario. This is 
a particularly attractive feature given the problems with 
fitting the details of a generic starforming model to high- 
and low-SFR galaxies. 

Why is it that the CO-RC-FIR relationship holds both 
globally and within galaxies to scales of a few hundred 
parsecs, but breaks down on smaller scales, scales of indi¬ 
vidual GMCs? Pressure mediation provides a natural scale 
for the relationship: the pressure scale height in galactic 



Fig. 8 . P for the Milky Way. 


disks. On larger scales, hydrostatic pressure is determined 
by the slow variation of the gravity provided by the stars, 
on size scales equal to or larger than the thickness of the 
gas layer in a galactic disk. In the vicinity of the Sun, this 
scale is about 250 pc, but the thickness varies with galac¬ 
tic radius. Superposed on this is a higher spatial frequency 
variation in the spiral arms; on these and smaller scales, 
local pressure effects can become important. 

Thus, our theory suggests that it is primarily HS pres¬ 
sure that drives the correlation by regulating both the 
fraction of gas in molecular form and the magnetic field 
strength. Secondary effects that can modify the basic cor¬ 
relation include (but are not limited to) spiral density 
waves, local star formation effects, and CR propogation. 

5.4. Spectral Index Dependence on Hydrostatic 
Pressure: A Leaky Box Model 

For our model to work requires a weak dependence of the 
CR electron density, Nq, on HS pressure (§ 5.3). The ordi¬ 
nary minimum energy argument instead accords equipar- 
tition to the energy of the relativisic electrons and that of 
the magnetic field, so that they both should scale with P. 
This argument would lead to a RC-CO correlation much 
steeper than observed: /rc cx whereas we instead ob¬ 
serve an exponent that is lower by an order of magnitude. 
Furthermore, it is commonly accepted that CR electrons 
diffuse from their injection sites. Thus, the CR electron 
density at a given location in the disk is the result of the 
balance between different processes such as the CR injec¬ 
tion rate, the electron confinement time and energy losses. 

To first approximation one might expect that the CR 
injection rate should be proportional to the SFR surface 
density, and thus, to the FIR emission. The FIR-CO cor¬ 
relation (Figure 121) would then lead to a direct scaling of 
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Nq with Ico- However, this in contrast with the hydro¬ 
static pressure model which predicts a weak dependence 
of No on Ico (see previous section). 

Let us consider the simple hypothesis that a galaxy, 
or any star-forming region it contains, behaves essentially 
like a leaky box which is continuously replenished by a 
constant injection of fresh particles with a power law en¬ 
ergy distribution at a rate dN{e)/dt = A(~^. There are 
two time-scales which are relevant for the CR electron 
density. One is the confinement time of the electrons in 
the star-forming region, tc, the other is the radiative life¬ 
time of the synchrotron electrons, tsyn- If tc ^ 4yn, OR 
electrons will leave the galaxy disk (or the star-forming 
region) before losing their energy. That is, we expect that 
as the hydrostatic pressure increases, the star-formation 
rate grows, and with it the FIR emission and the cosmic 
ray injection rate. However, the strong stellar winds asso¬ 
ciated with intense star-formation rates will facilitate the 
propagation of the “excess” cosmic rays from the disk to 
the halo. In this regime, the synchrotron emission spec¬ 
trum will be given by Eq. |21with No = Ate- Under the 
assumption of hydrostatic pressure regulation, we found 
that No is nearly independent of P (§ 5.3); thus we ex¬ 
pect tc oc 1/H, and we expect that the observed spectral 
index is about the same as the spectral index of the in¬ 
jected particles. 

As the hydrostatic pressure decreases, tc will progres¬ 
sively increase. When tc ^ tsyn, the CR electrons will lose 
all their energy inside the parent star-forming region and 
the synchrotron spectrum will be given by 

Inc oc NoB^^+^y^iy-^/^ , (6) 


with No = AjP"^. In this asymptotic regime, the spec¬ 
tral index of the synchrotron radiation emission is 0.5 
steeper with respect to the injection spectral index; i.e. 
Inc oc for an injection spectral index of ao = 0.8 

(see Fig. 13 . 

The electron synchrotron lifetime^ at a given frequency 
scales as tsyn oc . Thus, under the hypothesis 

of energy equipartition oc P, we have igyn oc 
In regions of strong star formation, where tc tsyn, we 
have seen that Inc oc NoP^'^ and Nq = Ate', therefore 
Inc oc AteP^'^ and 


tc IrcP-°-^ 
AP-3/4 • 


(7) 


If we can further relate the particle injection rate A to 
the star formation rate (and more tentatively, to the FIR 
luminosity) as A = TtSFR = Ifir, then 


JlL oc ^p-015. (8) 

^syn IpiR 

We have seen in § 4.4 that the FIR emission may have 
a shorter scale length than the RC and CO emission, at 

^ The characteristic radiative age for electron emitting at 
1 GHz in a 10 fiG magnetic field is ^ 3 x 10^ yr. 



Fig. 9. Spectral index between 327 MHz and 1.4 GHz as a 
function of tc/Gyn for a leaky-box model with an injection 
spectral index ao = 0.8, see text. 


least for galaxies with bright CO and relatively strong star 
formation. According to Eq. |3and Eig.|3 we expect that 
an increase of the RC/EIR ratio with radius will result in 
an average radial steepening of the synchrotron spectrum 
in galaxy disks. In Fig. cni we present the spectral index 
images of IC 342 and NGC 5194, measured between 327 
MHz and 1.4 GHz. The 327 MHz images are taken from 
the Westerbork Northern Sky Survey and their resolution, 
about 54", is well matched to that of the VLA D-array at 
1.4 GHz. At these two frequencies, the contribution from 
thermal emission is negligible. Figure ^1 shows that the 
spiral arms are characterized by a ~ 0.8. This indicates 
that these intense star-forming regions suffer a significant 
GR electron leakage. According with the predictions of 
the hydrostatic pressure model, in these high-pressure re¬ 
gions most probably tc Gyn- The inter arm regions, or in 
general the underlying disks, are characterized by steeper 
spectra. The radially averaged spectral index profile show 
a monotonic steepening of the synchrotron spectrum with 
radius up to a ~ 1.3 — 1.4. This spectral behaviour sug¬ 
gests that as P decreases, the confinement time of the 
radiating electrons in the disk became larger than Gyn- 
It is worth noting that Niklas & Beck (1997) used a 
similar argument to test the phenomenological model of 
Helou & Bicay (1993). Niklas & Beck concluded that on 
global scales, this model failed to produce the expected de¬ 
pendence on spectral index for galaxies with different ra¬ 
diation fields and magnetic field strengths. Indeed, the lin¬ 
earity of the global RG-FIR correlation, Sne oc S'fir indi¬ 
cates that (tc/Gyn) is, on average, constant. This explains 
why a relation between spectral index and magnetic field 
is not seen on global scales. However, we have extended 
this idea to local scales, within galaxy disks. According to 
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Figs. El and cni we expect tc to be some fraction of tgyn 
within disks. 

To summarize, we have used observationally deter¬ 
mined relationships between the molecular gas surface 
density and the hydrostatic pressure to express the CO 
and RC intensities in terms of P. Both quantities scale 
with pressure since high-pressure regions tend to be more 
molecular and, under the assumption of energy equipar- 
tition between the turbulent energy of the gas and the 
magnetic field energy, their synchrotron radiation ampli¬ 
fied. For the hydrostatic pressure model to work requires 
a weak dependence of the cosmic ray electron density on 
pressure in disk galaxies. The observed radial gradient of 
cosmic ray electrons in the Milky Way appears to sup¬ 
port this hypothesis. If cosmic ray leakage is taken into 
account, we show that this condition is compatible with a 
scenario in which the cosmic ray confinement time of the 
synchrotron electrons is inversely proportional to the SFR. 
The spatially resolved spectral index images of IC 342 and 
NGC 5194 are in general agreement with this prediction. 

6. Conclusions 

We have presented a comparison of the global RC, FIR, 
and CO in 24 BIMA SONC galaxies and a detailed study 
of the high-resolution (6") CO-RC correlation in a subset 
of 9 sources. Our main conclusions are as follows: 

(1) The global correlations between CO-RC, FIR-RC, 
and FIR-CO are comparably tight, and we cannot identify 
any of the relations as being a fundamental one. 

(2) At 6" resolution, for CO-bright galaxies, the CO- 
RC correlation is as tightly correlated as its global value, 
down to ~ 100 pc size scales. The correlation between the 
CO and RC emission is linear with a dispersion that is 
less than a factor of two. 

(3) Despite the low dispersion in the correlation, there 
are systematic variations in /RC • The organized struc¬ 
tures in gco/RC ^re found along spiral arms and on size 
scales much larger than the resolution of the observations. 

(4) Combining the above two points suggests that, 
down to 100 pc scales, there does not appear to be a char¬ 
acteristic size scale where the CO-RC correlation breaks 
down. 

(5) Unlike the FIR-RC correlation, there is no system¬ 
atic trend in CO-RC as a function of radius in galaxies. 
We speculate that high-resolution Spitzer observations will 
show detailed spatial correlation between the FIR and RC 
and CO. We also expect that Spitzer will show that the 
dominant variation in the FIR-RC correlation is across 
spiral arms, rather than in a radial trend. 

(6) We propose that the correlations between the RC 
and CO on large and small scales may be explained as 
a consequence of regulation by hydrostatic pressure in 
galaxies. We use observationally determined relationships 
between the molecular gas surface density and the hydro¬ 
static pressure to derive a relationship between the CO 
and RC intensities. Taken together with recent results 
linking the FIR emission with ISM pressure by Dopita et 


al. (2004), these scenarios link all three of the CO-RC-FIR 
emissions with the interstellar pressure, without invoking 
any explicit dependence on a star formation scenario. 

(7) Under the assumption of hydrostatic pressure reg¬ 
ulation, we use a leaky box model of cosmic rays to show 
that the spectral index between 327 MHz and 1.4 CHz 
depends on the ratio of the confinement time of the elec¬ 
trons in the star-forming region, tc-, to their synchrotron 
lifetime, tsyn- We expect that in regions of high hydro¬ 
static pressure, tc Cyn and a « 0.8. As the pressure 
decreases, the electrons decay before escaping the disk, 
tc ^ tsyn, and the spectral index steepens to a ~ 1.3. 
We show spectral index maps of IC 342 and NCC 5194 
that reflect this kind of dependence, with a « 0.8 in the 
spiral arms and a monotonic steepening with radius up to 
a w 1.3- 1.4. 

Acknowledgements. The authors thank an anonymous referee 
for the valuable comments and suggestions. We thank Hans 
Bloemen, Alberto Bolatto, Al Glassgold, Anne Green, Wim 
Hermsen, David Hollenbach, John Lugten, Chris McKee, Andy 
Strong, and Peter Teuben for helpful discussions. We are grate¬ 
ful to Rainer Beck, who provided us with his M51 radio¬ 
continuum image. Erik Rosolowsky kindly produced Figure |H| 
for us. T.T.H. is grateful to M. W. Regan for a helpful ex¬ 
change regarding the CO scale lengths. This work is partially 
supported by NSF grant 0228963 to the Radio Astronomy 
Laboratory at UC Berkeley. M.M. wish to thank the Berkeley 
Radio Astronomy Laboratory for their hospitality during the 
development of this work. 

Appendix A: Notes on Individual Galaxies 

IC 342. At d=3.9 Mpc, IC 342 is the nearest galaxy in 
the sample, and the high-resolution images provide linear 
resolutions of ^100 pc. (A recent Cepheid distance deter¬ 
mination by Saha, Claver, & Hoessel 2002 suggests that 
d=3.3 Mpc, which would provide slightly improved lin¬ 
ear resolution.) A prominent (though kinematically weak) 
north-south bar dominates the CO emission in the inner 
2 ' of this starburst source, with two spiral arms at larger 
distances north and south of the bar (Crosthwaite et al. 
2001). Within the central bar, the scatter of the corre¬ 
lation is higher than it is at larger distances. The most 
significant deviation from the correlation is found in the 
inner 15" where the peak of the RC emission is ~ 10 times 
higher with respect to the CO emission. This starburst 
galaxy may be a large-scale analog to starforming regions 
in the Milky Way, which do not exhibit the RC-FIR cor¬ 
relation on size scales of a few hundred pc (e.g. Boulanger 
& Perault 1988). However, it is relevant to note that on 
larger (55", or 1.0 kpc) scales, the CO-RC correlation re¬ 
mains tight. 

NGC 4258. The radio emission in this unusual source 
is associated with two, large-scale “anomalous arms” not 
seen in the optical but detected first in Ha (Courtes & 
Cruvellier 1961). In the inner galaxy, the radio emission 
is clearly associated with jets, presumably related to the 
ACN activity associated with the supermassive black hole; 
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Fig. 10. Spectral index images (left panels) and radial averaged profiles (right panel) for IC 342 and NGC 5194. The 
spectral index is measured between 327 MHz and 1.4 GHz. Only pixels with an intensity level above 3 ct at both 
frequencies have been considered. 


however, at larger distances, the nature of the radio emis¬ 
sion is not so clear (e.g. Krause & Lohr 2004 and references 
therein). A detailed discussion of this unusual galaxy is be¬ 
yond the scope of this paper, and a proper discussion of 
the RC-FIR-GO correlations in this source would warrant 
a separate study. However, we include NGG 4258 here and 
emphasize that, even in this highly unusual source, the 
GO-RG correlation is very good over those regions where 
both GO and RC are present. 

NGC 4414. The “low-resolution” RC emission from 
this source has been convolved from 14" resolution up to 
55", and at 55", the RC emission as well as the CO emis¬ 
sion are only marginally resolved. One must therefore be 
somewhat cautious about whether the 55" pixels shown 


in the point-by-point comparison in the upper right panel 
of Figure 1(c) are truly independent. However, at 6", the 
CO and RC structures of this galaxy are revealed as a 
ring with a diameter of ~ 5 kpc. NGC 4414 is a flocculent 
galaxy, which means that it lacks prominent spiral struc¬ 
ture in optical light, though Thornley (1996) found that 
K' imaging revealed low-level spiral structure. However, 
of four sources with newly-revealed spiral arms presented 
by Thornley, NGC 4414 remained “the most flocculent of 
the sample in K' emission”. The gco/RC image shown in 
Figureiniappears to be nearly featureless and shows a dis¬ 
persion of only 0.09; that is, the CO/RC intensity ratio 
varies by only a factor of 1.2. On the other hand, the dy- 
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namic range of emission in the CO and RC varies by only 
slightly more than one order of magnitude. 

NGC 4736. The molecular gas in this nearby source 
is distributed in a central molecular bar and in two tightly 
wound spiral arms which are associated with a bright ring 
of star formation (Wong & Blitz 2000). A background ra¬ 
dio continuum source, apparent in the upper left panel of 
Figure l{d), has been omitted from the analysis in this pa¬ 
per. As discussed in §4.3, there appears to be an enhance¬ 
ment in qco/RG (FigureEJ associated with the northwest¬ 
ern spiral arm; the enhancement may be a result of the 
starburst there. It also appears that qco/RC is larger at 
the bar ends, at roughly ±15" from the galaxy center, and 
along the “arclike” extensions from the bar that probably 
are associated with the beginnings of the spiral structure 
(Wong & Blitz 2000). 

NGC 5005. Two background continuum sources are 
apparent in the low-resolution RC images of NGC 5005; 
these have been omitted from all analysis. Like NGC 4414, 
the 55" observations show that the CO and RC is only 
slightly resolved. At 6" resolution, the central unresolved 
RC emission is associated with a nonstellar active galac¬ 
tic nucleus (AGN). Outside the nucleus, the agreement 
between the RC and CO emission is excellent, as can be 
seen both in Figure 1(e) as well as in FigureEl the gco/RC 
image is nearly featureless and has low dispersion. 

NGC 5033. Like NGC 5005, NGC 5033 has an un¬ 
resolved AGN at its nucleus, and its qco/RC image is 
nearly featureless. An enhancement in qco/RC is seen at 
the southern end of the CO emission; however, this CO 
emission is detected only at a ~ 3 (t level and might be 
spurious. Background continuum sources are apparent in 
both the high- and low-resolution data sets. 

NGC 5055. NGC 5055 is the prototypical flocculent 
spiral galaxy, which (like NGC 4414) was revealed to have 
regular spiral structure in AT'-band imaging (Thornley 
1996). At 6" resolution, the RC and the CO emissions 
appear nearly featureless, though an inspection of the CO 
image and especially the qco/RC image (Figure EJ show 
structures that follow very closely what appear to be inner 
spiral arms in Thornley’s AT'-band image. Nonetheless, the 
dispersion in 9co/RC Is only 0.11. 

NGC 5194. The Whirlpool galaxy M 51 is the pro¬ 
totypical grand design spiral galaxy, with two prominent 
spiral arms where the molecular gas is generally confined. 
The arms wind in towards the center of the galaxy but re¬ 
main in two distinct peaks; the molecular gas lacks emis¬ 
sion right at the nucleus. The RC at the nucleus is domi¬ 
nated by the AGN activity. The qco /RC image is remark¬ 
able in that it shows strong enhancements along the spiral 
arms, as discussed in §4.3. 

NGC 6946. This nearby starburst galaxy has an elon¬ 
gated, bright north-south molecular feature in the inner 
^ 50" and wispy arms at larger radii. The g'co/RC image 
(Figure EJ and profile (Figure li) show an enhanced ra¬ 
tio within the central region, apart from the compact HII 
region nucleus, which shows relatively enhanced RG emis¬ 


sion. Like IC 342, it may be that this region is a large-scale 
analog to active starforming regions in the Milky Way. 
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